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Abstract 

A new kinetic approach is developed and a quasi-practical gas is de- 
fined to which the new approach can be applied. One of the advantages 
of this new approach over the standard one is direct calculability in terms 
of today's computational means. 

PACS 51.10.+y - Kinetic and transport theory of gases 



1 Introduction 

The impressive progress of nonlinear science in the last several decades revealed 
a large number of structures and mechanisms that were novel and fascinating 
to almost all of us. This has, in a sense, suggested that several aspects of our 
existing statistical theory, whose first kinetic equation was proposed more than 
a century ago, need certain kinds of 'upgrade' or 'renewal'. 

In some of our relatively recent papers [1, 2], we were tempted to analyze 
the standard kinetic theory from different perspectives. It was argued that (i) 
in various artificial and realistic situations, distribution functions of gases may 
have very unconventional structures: they can be, for instance, discontinuous at 
every and each spatial point and thus the usual differential apparatus becomes 
inapplicable (after all, there are indeed many structures and mechanisms in 
nature to which the use of the usual differential apparatus is quite limited) ; and 
that (ii) while the left side of the Boltzmann equation is symmetric in terms of 
partial derivatives with respect to the position and velocity, the right side of the 
equation is constructed entirely in the velocity space (the position coordinates 
merely serve as inactive parameters), which is incomprehensible especially in 
the mathematical sense. 

The major objective of this work is, however, rather practical: to define 
a kinetic gas and, at the same time, to find out an algorithm which can be 
used to calculate the behavior of the gas rather completely. It is hoped that 
if such study gets established somehow, further studies, either concerning the 
foundation of the existing statistical theory or concerning the development of 
more general treatment, will be inspired and promoted. 

The structure of this paper is the following. In Sect. 2, we propose our 
working model, a gas leaking out of a container, resembling the cavity model 
for the black-body radiation. In Sect. 3, the zeroth-order distribution function 
of the leaking gas is formulated and, in passing, it is shown that distribution 
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functions of real gases may have radically discontinuous structures untreatable 
for the standard theory. In Sect. 4, the collisional correction is investigated with 
help of a methodology that is slightly different in form but much different in 
concept from the standard treatment. In Sect. 5, we further comment on why 
the distribution function averaged over finite velocity solid-angle ranges needs 
to be introduced. Sect. 6 summarizes the paper. 

2 A leaking gas 

At the kinetic level, one of the customary conceptions about calculating a prac- 
tical gas is to solve the Boltzmann differential-integral equation with help of a 
difference scheme. Known difficulties associated with this conception may be 
summarized as follows. 

1. There are seven variables: time, three spatial coordinates and three veloc- 
ity components. If we divide each variable into N intervals, 'the degrees 
of freedom' of the system are iV 7 . To reveal the true properties of a 
nonequilibrium phenomenon, such as those related to a turbulence, N has 
to be rather large and N 7 has to be terribly huge, so that no today's 
computational means really helps. 

2. The complex nature of the collisional operator in Boltzmann's equation 
makes the situation worse. 

3. Boundary conditions and initial conditions impose extra problems. As to a 
differential equation, well specified boundary and initial conditions usually 
mean that a uniquely and clearly defined solution can be constructed, say, 
from a difference scheme. (The uniqueness may not be truly desired in 
view of the fact that bifurcations take place in nature.) As to an integral 
equation, boundary and initial conditions may be specified rather loosely 
and it leaves us certain room to 'manipulate'. Boltzmann's equation is a 
differential-integral equation and it is not entirely clear in which direction 
we should go. 

Due to these difficulties (and possibly many more), almost no realistic problems 
have been fully treated, let alone conclusive comparison between the kinetic 
theory and realistic experiments. 

Interestingly enough, the aforementioned difficulties, at least some of them, 
are not intrinsic to the dynamics that we wish to study. In some sense, if the 
standard approach had not dominated our mind too strongly, we might have 
already had workable alternative approaches, at least for some special cases. In 
what follows, we shall try to substantiate this viewpoint. 

We consider a gas consisting of hard spheres and confined to a closed con- 
tainer. The interactions between the walls and particles, and between the par- 
ticles themselves, make the gas solidly and quickly in equilibrium. Namely, the 
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probability of finding particles in a phase volume element drdw — dxdydzdv x dv y dv z 
takes its Maxwell form 

fM 55 n ° 2 exp ("^) ' (1) 

where m is the mass of particle, k the Boltzmann constant, T the temperature 
of the gas and n the particle density in the container. 



Figure 1: Schematic of a leaking gas. 

We then suppose that there is a small hole on the wall of the container, 
as shown in Fig. 1. The aim of our formalism is to determine the distribution 
function of the leaking gas. For convenience of discussion, it is further assumed 
that (i) the hole is indeed small so that the leaking is relatively slow and the 
relevant distribution function can be considered to be time-independent; (ii) 
since the regions far from the hole are kept in the vacuum state (by a pump 
for instance), no incoming particles need to be taken into account; and (iii) the 
leaking gas outside the container is so dilute that the zeroth-order solution can 
be determined by the collisionless trajectories of particles while the first-order 
correction can be formulated by assuming the particles of the leaking gas to 
collide with each other once and only once. 

3 The zeroth-order distribution function 

If we ignore the particle-to-particle interaction and ignore the gravitational force 
acting on each particle, the collisionless motion of every particle of the leaking 
gas is conceptually simple — moving along a straight line. However, it is still 
necessary and interesting to express the corresponding distribution function in 
a proper mathematical form. 

We start our discussion with the conventional theory. The theory states that 
the distribution function satisfies, with collisions ignored completely, 



where F represents the external force. Equation (2) is sometimes termed the col- 
lisionlcss Boltzmann equation [3]. A natural notion related to such terminology is 
that any procedure, a difference scheme for instance, if applicable in solving the 
regular Boltzmann equation, must be applicable in solving this reduced equa- 
tion. Strangely, complicated and subtle issues arise from this 'natural' notion; 
and we shall discuss these issues at the end of this section. 

To formulate the zeroth-order distribution function, the following 'slightly 
different' approach is helpful. Rewriting Equation (2) along a particle path in 
the phase space, we arrive at the path invariance[3, 4] 



As far as our leaking gas is concerned, this means /| pa th = /m, where f M has 
been defined by expression (1). 



Figure 2: The velocity distribution at the position r. 

Then, it seems that the zeroth-order distribution function on the right side 
of the hole is uniformly identical since any spatial point there is reachable along 
a path radiating from the inside of the container (the external force is zero 
and the paths in the phase space and in the spatial space are the same). This 
notion is, however, rather misleading. By moving together with a particle of 
the leaking gas, an observer easily realizes that the true particle density around 
him becomes lower and lower. With this realization in mind, it can easily be 
found that the distribution function of the leaking gas is strongly limited in the 
velocity space: the farther from the hole the stronger the limitation is. Namely, 
as shown in Fig. 2, the true implication of (3) is 



where Ail is a solid-angle range whose size is defined by the solid-angle range 





Jm the direction of v within ACl 
otherwise 



(4) 



AS cos 9 



with r = |r 



(5) 
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and whose direction is defined by the position vector r (with the origin at the 
hole). Heuristically and intuitively, it can be said that the velocity distribution 
at every point there is a sting-like function. In almost all the regions on the 
right side of the hole, where r is relatively large or cos 9 is relatively small, the 
stings tend to be infinitely acute. 

Though (4) and (5) correctly describe the distribution function, there are 
some kinds of inconvenience in applying them in analytical calculations. For 
practical and theoretical reasons, we wish to express the distribution function 
/(r,v) directly and explicitly in terms of r and v. 

To get such expression, we adopt the picture that the sting-like velocity 
distributions of the leaking gas are indeed infinitely acute, either by assuming 
the hole to be truly small or by assuming the interested region to be truly distant. 
Under this understanding, we can use a 5-function to reflect the limitation of 
the distribution function in the velocity space and replace (4) and (5) by 



where f2 v is the solid angle in the direction of v and fl r is the solid angle in 
the direction of r. Expression (6) enables us to do analytical calculation easily. 
For instance, with help of it we may compute the zeroth-order particle density 
outside the container by 



which is quantitatively consistent with our intuitive notion about the particle 
density outside the container. 

Before leaving this subject, we go back to discuss several usual concepts 
concerning the collisionlcss and regular Boltzmann equations. It is commonly 
believed that the collisionless Boltzmann equation (2) is completely equivalent 
to the path-invariance expressed by (3). In our view, this equivalence is just a 
formal one. For one thing, expression (3) makes sense strictly along a particle 
path, which is governed by the collisionlcss motion equations of single particle; 
while Boltzmann's equations are supposed to be solved by a difference scheme 
with help of boundary and initial conditions in which trajectories of individual 
particle do not play any role. For another, unconventional topology and dynam- 
ics are related to (3). Distribution functions described by it, such as (4) and 
(6), can be shaped like stings; more than that, these stings 'spread' out along 
paths of particles, become sharper and sharper when spreading, and continue 
to spread out after being infinitely sharp. (Distributions produced by bound- 
aries can be even 'sharper', see Sect. 5.) Boltzmann's equations, being ones 
containing differential terms, are not compatible with these 'radical' things. 



/( r , v ) = g(r,v)S (O v - n r ) 



(6) 



with 
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4 The collisional correction 



In this section, we try to deal with collisions. To make our discussion involve less 
details, we shall only formulate the first-order correction. Namely, it is assumed 
that the particles expressed by (6) will collide once and only once (although 
further extension can be made with no principal difficulty). 




Figure 3: A particle detector is placed in the leaking gas. 

What can be measured in an experiment should be of our first concern. For 
this reason, we consider a particle detector placed somewhere in the leaking gas 
and try to determine how many particles will be recorded by the detector, as 
illustrated in Fig. 3. There are several essential things worth mentioning about 
the detector. 

1. The opening allowing particles to enter the detector is considered to be 
sufficiently small. Or, in the mathematical language, we regard the area 
of the opening, denoted by AS, as an infinitesimal quantity throughout 
this section. 

2. Without specifying concrete structure and mechanism of the detector, it is 
assumed that every particle that enters the detector and is in the velocity 
range 

Av = Aw 2 Afl, (8) 

where Av and Af2 are predetermined before the measurement, will be 
recorded, and every other particle will not. 

3. While AS and Av are allowed to shrink to zero, we define AQ as a finite 
solid-angle range. As will be seen, the discrimination against AO is taken 
almost entirely by necessity. 

If the detector works as described above, do we know the distribution func- 
tion at the entry of the detector? The answer to it is almost a positive one. If 
AN is the number that the detector counts during At, the particle density in 
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the phase volume element (ASvAt)(Avv 2 AH) is 

AN 

The form of (9) illustrates one of the most distinctive features of this approach: 
it tries to calculate the distribution function directly rather than to formulate a 
dynamic equation concerning how many particles enter and leave a phase volume 
element. Apart from other advantages, this brings to us a lot of convenience in 
computational terms. 

Noticing that only the collisions taking place within the shaded spatial cone 
— Afl, which is equal and opposite to the velocity cone Ail, can possibly con- 
tribute to AN, as shown in Fig. 3, we call the region inside —Ail the effective 
cone. Since this effective cone is a finite one (Ail is finite), our primary task is 
to divide it into many subregions, denoted by dr' (the origin of the coordinates 
is still at the center of the hole as in Sect. 3), and to calculate how collisions in 
each of the subregions give contribution to AN. 

Consider that two beams of identical particles (but still distinguishable in 
terms of classical mechanics) 

f(t'y,v')dV and /(tVi.rOdvi (10) 

collide within dr' and at time t', producing particles with velocities v and vi 
respectively. It is noted that there is a time delay 

t - t 7 = l r ~ r I with v = \v\ (11) 

v 

between the time of collision and the time of detection. Since the leaking is 
assumed to be relatively slow (or the gas inside the container is supplied by an 
external gas source), we think of our problem as a time-independent one and 
pay no more attention to the time variable hereafter. To better observe the 
collisions, we define 

{ V ! + V hf, and { v + v 1 = 2c 

[ v — Vj = 2u [ v — vi = 2u v ' 

That is to say, by virtue of the conservation laws of energy and momentum, 
c = c' represents the velocity of the center-of-mass and u = |u| = |u'| stands for 
the particle speed relative to the center-of-mass. Examining the beam-to-beam 
collision in the center-of-mass frame, we find that the differential number of the 
colliding particles is 

[dr'f(v', v')dv'} [/(r', vi)dvi] [2tKr c (u', u)dil c At] , (13) 

where il c is the solid angle between u' and u and a c is the cross section associated 
with particles emerging within the solid-angle range dfl c . By integrating (13) 
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and getting help from (6), the right side of (9) is equal to 



J —An JAvAUo Jo Jo (|r — r'\ 2 AQ, v)(v 2 AvAQ,) 

where Af2 is the solid-angle range formed by the point dv' (as the apex) and the 
detector opening AS (as the base), and the subindex AiAf2 there states that 
only particles that can be recorded by the detector will be taken into account. 
Since the size of AOo is 'much smaller' than that of AS1, every particle starting 
its journey from the effective cone and entering the detector will be treated as 
one within AO,. With help of the variable transformation from (v', v[) to (c', v!) 
and finally to (c, u), we rewrite expression (14) as 

f dr'f dn c fdc r d J^-\\J\\- 9 (c + u) 9 (c-u) t 
J— An JAvAQo Jo J-c (\r -r'\ 2 An v){v 2 AvAn) 

in which the Jacobian between the variable transformation is 

||J|| = f^=2. (16) 
o(c, u) 

In view of the symmetry of the cross section there, we have 

J du--- = 2 J du---. (17) 



Av '•' ' 

Figure 4: The relation between the velocity element v 2 AvAflo 
and the velocity element u 2 dud£t c . 

And, by investigating the situation in the velocity space shown in Fig. 4, the 
following relation can be found out: 

/ u 2 dudn c ■ ■ ■ « v 2 AvAfl ■■■■ (18) 

J AvAQo 

Therefore, the first-order distribution function at the entry of the detector is 



/(r,v,AQ) = -L- / dv' T dc 



8cr c (u', u)g(r', c + u)g(r', c-u) 



(19) 
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in which, although u = \c — v|, u < c needs to be ensured for the integral to 
make sense and the direction of u' is the same as that of c' or c. 

Since we have, at the beginning of this section, assumed Af2 to be a finite 
solid-angle range, expression (19) is nothing but the distribution function av- 
eraged over Af2, which is, in a sense, still different from 'the true and exact 
distribution function' there. In the next section, we shall come back to this is- 
sue and find out that the true and exact distribution function is actually beyond 
our reach. 

If interested, readers may test the resultant formalism with aid of realistic 
and computational experiments or develop it to cover more practical and more 
complicated cases. This is one of the major purposes of the present paper. 

5 Discussion 

In this section, we wish to further justify one of our introduced concepts, the 
distribution function averaged over finite solid-angle ranges of velocity. 




Figure 5: Particles refleted by a surface. 

Let's look at the path-invariance theorem expressed by (3) again. Although 
the path-invariance is not equivalent to the collisionlcss Boltzmann equation 
(see the remark at the end of Sect. 3), the theorem can still be thought of as a 
part of the conventional theory since it is formally the same as the collisionless 
Boltzmann equation and it is consistent with the usual notion that distribution 
functions of real gases are continuous or, at least, piecewise continuous. Then, a 
crucial question arises. Arc all distribution functions of realistic gases continuous 
or, at least, piecewise continuous? If this is the case, we will be able, at least 
in principle, to determine the exact distribution function. To see what the 
real situation is, we examine the case illustrated in Fig. 5, where particles are 
moving toward a surface with a definite velocity (or, equivalently, let the surface 
move toward the particles). Since the surface is not uniform geometrically and 
physically, we need to, in order to formulate the reflected particles, divide it 
into many small surface elements, denoted by AiS herein. It is easy to see that, 
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corresponding to one of AiS the reflected particles are like ones emitted from a 
point particle source and the distribution function, at a position r around the 
surface, is 

A»/= f \ V)AiS . 2 S(Q v -n r . Fi ), (20) 
47r|r — Y t \ z 

where rj represents the position vector of AiS and 77* (v) is the reflection (emis- 
sion) function that can be determined by experiments [5] . As has been illustrated 
in Sect. 3, the distribution function expressed by (20) is an infinitely thin 'sting' 
at every spatial point, and the particle density along a particle path is no longer 
invariant. The total distribution function at r is 

/ = £^/ = £S^^v-o r - ri) . (2D 

i i 1 

It is then found that, if the true and exact distribution function is indeed of 
our concern, the <5-functions in (21) will stay there and no regular functions 
can be used to replace them because each of f2 r _ r . in it points in a distinctive 
direction, which means, in a more vivid language, the velocity distribution at 
every spatial point is shaped like an infinite number of infinitely thin stings 
(similar to functions dealt with in the studies of fractals). Without introducing 
the distribution function averaged over velocity solid-angle ranges, we will have 
great difficulty in setting up a calculable theory. 



AO 




(a) (b) 
Figure 6: The effective region above the detector: (a) If AO is 
fixed and AS* shrinks to zero, it is cone-shaped, (b) If AS is fixed 
and Af2 shrinks to zero, it is cylinder-shaped. 

Treatment of particle-to-particle interactions is always an indispensable chap- 
ter of statistical mechanics. While referring readers to elsewhere[l, 2] for analy- 
ses of perplexing problems in the textbook methodology, we shall here comment 
on our own treatment of collision presented in the last section. At first glance, 
if Afl there were infinitely small, both the denominator and numerator of (19) 
would tend to zero simultaneously and the total expression would remain finite 
and valid. A careful inspection, however, tells us that there are two 'competing' 
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quantities Aft and AS. If we fix Aft and let AS approach zero, the effec- 
tive region is a cone-shaped one, shown in Fig. 6a, and the related distribution 
function is, as has been formulated, the one averaged over the finite solid-angle 
range Aft. Whereas, if we fix AS and let Aft approach zero, the situation is 
rather different: the effective region becomes a cylinder-shaped one, shown in 
Fig. 6b, and the related distribution function will be, if formulated, the one 
averaged over the finite spatial area AS*. In both the above situations, each of 
the results is related to an integration, each of the integrations is taken over 
a region infinitely-extended in space. (Even if more collisions are taken into 
account the above assertion still holds its significance since particles from the 
remote regions can still come freely in terms of probabilities.) There is no way 
to ensure the continuity of each formalism let alone the consistency between the 
two formalisms. 

The discussions in this section clearly show that, in dealing with both col- 
lisionlcss dynamics and collisional dynamics, a finite range of certain variable 
(whether in the position space or in the velocity space) and the average over this 
range have to be employed. If the approach is inherently and completely that of 
differentiation, which means that the exact distribution function is of concern 
and every position and velocity range must shrink to zero, insurmountable con- 
flicts surface automatically in one way or another. The repetitive appearances 
of conflict remind us of the situation in quantum mechanics where the accu- 
rate position and momentum cannot be determined simultaneously. Possible 
implication of this issue has yet to be explored. 

6 Summary 

We have proposed a gas as our working model and formulated a feasible method 
to calculate it. The formalism shows that real gases, at least certain types of 
them, can be calculated at the kinetic level. It is expected that such calculations 
will soon be compared to realistic or computational experiments. 

In this paper, collisional effects were investigated under the assumption that 
each particle involves, at most, one collision. If more collisions are taken into 
account, the trajectory of a particle will be very much like that of the Brownian 
motion or, in another sense, similar to that appearing in the logistic map of 
nonlinear studies. Although the extensions in this direction have been tried[6, 7] , 
it is not appropriate to comment on them before this work is admitted by more 
in this community. 

Helpful discussions with professors Hanying Guo, Ke Wu and Qiang Chen 
are greatly appreciated. This paper is supported by School of Science, BUAA, 
PRC. 
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